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We propose an experiment in which an entangled pair of optical pulses are propagated through 
non-uniform gravitational fields. A field operator calculation of this situation predicts decoherence 
of the optical entanglement under experimentally realistic conditions. 
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The effort to try to find a consistent way to combine 
quantum theory and gravity is made particularly difficult 
by a lack of experimental indicators. Typically, situations 
in which competing approaches make testable predictions 
involve experimental scenarios far beyond the reach of 
current technology, such as interactions with black-holes 
[l|, worm- holes 0, highly accelerated frames Q, exotic 
gravitational potentials [4j or on extremely short length 
scales Q . Because the non-local correlations of entangled 
quantum systems contrast the local nature of relativity, 
a fruitful area to look for experimental indicators would 
seem to be the behavior of entanglement in non-inertial 
frames [f|. 

Here we consider bipartite entanglement in which the 
members of the entangled pair experience different grav- 
itational holds before being measured. In particular we 
examine the effect on optical entanglement of propaga- 
tion through varying gravitational fields. Using a novel 
approach we predict an apparent decoherence of the en- 
tanglement under experimentally realistic conditions. 

Our proposed set-up is shown in Fig.l. We consider 
time-energy entangled photons p} that are sent on two 
different paths to a common end-point. We contrast this 
with classically correlated pulses sent on the same paths. 
The photons begin and end their journeys in regions of 
space with the same gravitational field. On the journey 
one of the photons passes through a region of significantly 
different gravitational field. We suppress the transverse 
dimensions and consider only time, t and space in the 
longitudinal direction, x. We define a localized spatio- 
temporal mode with average frequency lj and wave num- 
ber k, centred on the space time point x,t via the mode 
operator: 



a(t, x) 



dt'dx'e^'-^Git' -t,x'- x)a t >, x ' (I) 



where at\ X ' ar e single space-time bosonic field an- 
nihilation operators with the non-zero commutators 

[at',x',4",x"] = S ( t ' - t")6(x' ~ %")■ G(t',x') is a nor- 
malised temporal and spatial wave function, defining the 
measurement and interaction bandwidths of the mode. 
Single photon states of this mode are then given by: 



\l) t , x = tf(t,x)\0) 



(2) 



where |0) is the electro-magnetic vacuum state. 

The time/ position representation is related by Fourier 
transform to the frequency and wave-number spectra of 




FIG. 1: Proposed experiment. Vacuum inputs at A are pro- 
cessed by a source S which produces two output modes. Mode 
1 propagates orthogonal to the gravitational field direction till 
it meets mirror ml and is reflected back to the detector at B. 
Mode 2 propagates anti-parallel to the field till it meets mir- 
ror m2 and it is reflected back parallel to the field to the other 
detector at B. The photo-currents from the detectors are com- 
bined on a multiplier to read out the coincidence current. 



the mode thus an alternative representation of Eq[T] is 



8]: 



a(£, x) 



dw'dk'e t ^' t - k '^G(uj' -u,k'- k)a u 



where the bar indicates Fourier transform of the relevant 
function and the a u i t i-' are single wave-number/frequency 
bosonic annihilation operators with the non-zero commu- 
tators [a u i,k', a u " ,k"\ = — v")8{k' - k"). 

Free evolution over a time interval s in a region of 
space with negligible gravitational field is generated by 
the unitary operator U s i which displaces the state in time 
by the amount s and in space by / = s, where we work 
in geometric units for which the speed of light is unity. 
Explicitly, the unitary can be written 



U, 



Sji = exp[i 



dtu"dk"(uj"s - k"l)al„ k „a^^] (4) 



ft is straightforward to confirm that U^ sl a{t, x)U St i — 
a(t + s, x + I) as expected for Heisenberg evolution of 
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the mode operator. We wish to generalize this expres- 
sion to treat radial propagation close to a massive body 
of mass M. For simplicity we will assume the body is 
non-spinning. The Schwarzschild metric in the radial di- 
rection for such a body is given by 



, n , 2M, , „ 

dr 2 =1 )dt 2 

r 



dr 2 



(5) 



where t is the time interval measured by clocks in a dis- 
tant inertial frame and r is the reduced circumference. 
However, noticing that the local proper intervals at a 
stationary "shell" frame at radius r are given by [9j 



ds — 
dl = 



2M , 
1 dt 

r 



dr 



2M 



(6) 



we can rewrite EqJS] as dr 2 = ds 2 — dl 2 . For free opti- 
cal propagation dr = 0, hence ds = dl and local evo- 
lution of the mode in a shell frame is still governed 
by Eqf5] provided we make the substitutions s — > ds 
and I — > dl. In addition, lj and k in the definition of 
the mode operator, EqJ31 should be interpreted as the 
locally measured shell values of the average frequency 
and wave-number. The global radial evolution can be 
evaluated by integration over the local evolutions such 
that the Heisenberg evolution becomes d(t,x,u),k) — > 
d(t + j ds, x + J dl,u>', k'), where lj' = g u> and k' — g^ 1 k 

with g = y/l - ^/v/l - Tf and r e an d r f the initial 
and final radii respectively. In the following calculations 
the detectors and source will be in the same shell frame 
such that rf = r e and the initial and final mode opera- 
tors will have the same values of the average frequency 
and wave-number. We now use these tools to solve the 
problem posed by Fig.l. 

We will refer to the shell frame in which the source 
and detectors are located as the SD-shell. We assume 
that photon mode 1 remains approximately in the SD- 
shell throughout its propagation, i.e. r = r e is a constant. 
The evolution of photon mode 2 is evaluated in a succes- 
sion of shell frames as the pulse climbs out of and is then 
reflected back into the gravitational well. The mirrors 
ml and m2 are assumed ideal reflectors across a band- 
width much larger than that of the photon modes. We 
set two conditions on the evolution: (i) the time intervals 
for the two modes to propagate to their respective mir- 
rors and back as measured in the SD-shell are equal and; 
(ii) the total time between the relevant initial states and 
the final detection events, as determined by a sequence of 
shell-frame observers along the light pulse paths, arc also 
equal. Condition (i) constrains the position of the mir- 
rors such that classical pulses launched simultaneously 
from the source will strike the detectors at the same time. 
Condition (ii) ensures equivalent unitary evolution along 
all paths. 



We use the Heisenberg Picture to evaluate the expected 
photon counting correlations by evolving the mode oper- 
ators from the photon counters at B back to the initial 
vacuum states at A and expressing them as functions of 
the input operators. In general we will get solutions 



a m ,i = fi(ai(s,l),a 2 (s / > l')) 
a m ,2 = hWX),a 2 {s"',l"")) 



(7) 



where the functions /, describe the mixing of the mode 
operators and the classical pump field by the source. The 
different time-space labels on the operators can arise due 
to the different evolutions produced by the gravitational 
field. We assume all these times precede the source in- 
teraction, i.e. s,s',s"s"' « t , where to is the time at 
which the interaction with S took place. As a result all 
the input operators act on the vacuum state. The rate of 
coincident detection events, C, is given by the expecta- 
tion value of the product of the photon number operators, 



of the two modes. Thus 



C = (nifia) = (00|aJi 



m : l L 



,l& m 2 dm,2 



|00) 



(8) 



In order to solve Eq[5]we must insert the correct source 
mode transformations and evaluate the relevant time- 
space labels for the mode operators given the conditions 
(i) and (ii). First we evaluate the integrated shell- frame 
time (1ST), a c (see Fig.l), of mode two as it climbs ra- 
dially out of the field, by integrating over a succession of 
shell frames as it climbs. We obtain 



dr 



1 - 



2M 



(9) 



where r e + h is the reduced circumference of mirror ni2. 
By symmetry the 1ST of the mode as it propagates back 
from the mirror to r e is also equal to a c . We also require 
the value of this same interval but as measured by clocks 



(SB) 



in the SD-shell, a. 
ds = dl we have the relationship 



From Eqs [6] and the condition 



dt 



dr 



Hence we obtain 



r (SD) 



2M 



dr 



(1- ™ 

r (SD) 



(10) 



(11) 



From condition (i) we require that a j — a c where a j 
is the time taken by mode one in the SD-shell (see Fig.l). 
Condition (ii) requires that o\ = 02 where 



a 1 = 2a s + 2a c + a vl 
a 2 = 2a s + 2a f + a v2 



(12) 



where the various other time intervals are defined in 
Fig.l. Substituting from Eqs [91 and [TT1 we find 



A 



— &v2 = 2(a f — a c 
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2M 

2(\/!-— >< 



rf?" 



1 _ 2M 



(13) 



and so we can write the space time labels of mode oper- 
ators evolved back via path one as 



Si — ll — to — CFvl 

and that for path two as 

s 2 = h = t - a v i - A 



(14) 



(15) 



For weak fields, 2M/r << 1, and heights much smaller 
than earth radius, h/r e << 1, we obtain the approximate 
expression 



A 



-h 2 M 



(16) 



Now we can solve for particular sources in Fig.l. Con- 
sider first that the source S produces classically corre- 
lated coherent pulses. The displacements introduced by 
the classical source are assumed matched to the quantum 
modes. Explicitly, the input/output relationship for such 
a coherent source applied to mode 1 is given by 



flout, l 



J f dt'dx'G(t' — t ,x' — x )(ai t t',x' 

+ G(t' - t ,x' - x )a) 
a ini x + a (17) 



where a is the coherent amplitude of a classical dis- 
placement field, and normalization of the spatio-temporal 
mode function has been used in going from the first to 
second line. Similarly for mode 2 



flout, 2 



a™, 2 + a 



(18) 



where condition (i) ensures that both modes are matched 
to the same classical displacement. Using Eqs [TJ] - [T5] 
we obtain the following expression for the evolved mode 
operators: 



a m ,i = ai(si,h) + a 
a m ,2 = a 2 (s 2 ,l 2 ) + a 



(19) 



The mode operators appear with different time space la- 
bels, but commute for all values, and act upon their re- 
spective vacuums to produce the null state. Hence using 
EqJ5]we obtain 



C=\a\ 



(20) 



As expected the classical correlation is unaffected by 
propagation through the varying gravitational fields. In 
general this result will hold for all separable correlations. 



Now we consider the source S in Fig.l to be entan- 
gling. In particular we consider the production of time 
energy entanglement from the vacuum inputs via para- 
metric down conversion. The input/output relationship 
for the down conversion is given by 



Clout, 1 
flout, 2 



flm.l 
«m,2 



(21) 



where \ is proportional to the classical pump field of the 
down converter, which again is assumed mode matched 
to the detection modes, and we assume x << 1. In this 
limit the interaction described by Eqs [21] correspond to 
the Schrodinger evolution |0,0) — > |0, 0) +x|l)t,x,i|l}i,x,2- 
This state is weakly entangled. A signature of the en- 
tanglement is that, although photon events are rare and 
random, they are perfectly correlated between the beams, 
i.e. the photons always arrive in pairs. Using Eqs [T4l[T5l 
and [21] we obtain the following expression for the evolved 
mode operators: 



a m ,i = ai(si, h) + xal(si,h) 
a m ,2 = a 2 (s2,Z 2 ) +Xa\(s2,l2) 



(22) 



Now the commutation relations between the mode oper- 
ators are changed in a non-trivial way. To second order 
in x w e obtain for the correlation 



C 



The commutator 



\x[ai(si,h),a\(s 2 ,l2 



(23) 



a 1 (s 1 ,l 1 ),a\(s 2 j2)) 

dt'dx'G{t' - s u x' - h)G{t' - s 2 ,x' - l 2 ) 



(24) 

is a real number between and 1 that depends explicitly 
on the mode functions. However, in general it will equal 1 
when A = and will equal when A >> 0. Thus in the 
presence of a uniform gravitational field we observe per- 
fect photon number correlations whilst in the presence of 
a non-uniform gravitational field effect of sufficient mag- 
nitude those correlations completely disappear. The en- 
tanglement appears decohered. 

To estimate the size of this effect we consider an exper- 
imental apparatus at sea-level where mode 2 is sent ver- 
tically a distance h before being reflected back to earth. 
Assuming the mode function is a Gaussian of the form, 



G(t-t',x-x') 



1 71 



ndtd x 

we obtain from Eqs [23] and l24l 



C = X "e d t <'. 



(25) 



(26) 



We estimate the intrinsic temporal and spatial uncer- 
tainties of a silicon photon counter to be 300 fs and 1 
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FIG. 2: Plot of the normalized coincidence rate, Cjv = C/x 2 , 
versus the height in metres, h, of the mirror m2 (see Fig.l). 
We use the parameters (all expressed in units of length for 
which c = 1) r e = 6.38 x 10 6 m, M = 4.432 x 10" 3 m, dt = 
10~ 5 m and d x = 10" 3 m. 



mm resectively. Hence, we set the standard deviations 
in Eq(2S]to d t — 10~ 5 to and d x = 10~ 3 m. These num- 
bers also roughly correspond to the pulse duration and 
pump spot radius for typical, single pass, pulsed para- 
metric down-conversion. The results of these numerical 
calculations are shown in Fig. 2. We find strong deco- 
herence (i.e. disappearance of coincidences) for h > 400 
kilometers. 



The novelty of the predicted effect should not be under- 
estimated. First note that although, because of the loss 
of photon correlations, we refer to this effect as decoher- 
ence, in fact the effect is in principle reversible by resend- 
ing (before detection) mode 1 along mode 2's path and 
vice versa. Secondly we note that energy conservation, as 
usually expected from down conversion, is only satisfied 
on average, not shot by shot. On the other hand in the 
current model both the down conversion pump and the 
gravitational field itself are treated as un-depleted clas- 
sical energy reservoirs, so a definitive statement on this 
issue is beyond the scope of the present model. Thirdly 
we anticipate that more unusual evolutions may arise 
for strongly entangled qubit states as suggested in Ref 
[l(J. Treatment of such situations with the same rigour 
as used here would require consideration of highly non- 
linear Heisenberg evolutions that are, again, beyond the 
scope of the present calculations. 

We have studied the effect on optical entanglement of 
evolution through varying gravitational fields. We have 
predicted a decoherence effect that should be observ- 
able under experimentally achievable conditions. We be- 
lieve the outcome of such an experimental investigation 
could have considerable implications for the unification 
of quantum mechanics and general relativity. Further- 
more, the predicted effect, if observed, would represent a 
new phenomenon with major consequences for quantum 
physics in general and quantum information in particu- 
lar. 
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